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Color- Coulomb Force Calculated from Lattice Coulomb Hamiltonian 
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The static color-Coulomb potential is calculated as the solution of a non-linear integral equation. This equation 
has been derived recently as a self-consistency condition which arises in the Coulomb Hamiltonian formulation 
of lattice gauge theory when the restriction to the interior of the Gribov horizon is implemented. The potential 
obtained is in qualitative agreement with expectations, being Coulombic with logarithmic corrections at short 
range and confining at long range. The values obtained for the string tension and Ajjg are in semi-quantitative 
agreement with lattice Monte Carlo and phenomenological determinations. 



1. INTRODUCTION 

The lattice Coulomb-gauge Hamiltonian has 
recently been derived from the transfer matrix of 
Wilson's Euclidean lattice gauge theory jl). The 
physical configuration space (no Gribov copies) is 
restricted to the fundamental modular region A 
of the minimal Coulomb gauge — i.e. to the set 
of absolute minima of a Morse function on the 
gauge orbits. This restriction is implemented by 
the effective Hamiltonian 



H 



eff 



K 



coul 



7o G 



(1) 



where H CO ui is a lattice analog of the Christ-Lee 
Hamiltonian and G is the (3-dimensional) 
horizon function Here 70 is a new thermo- 
dynamic parameter whose value is determined 
by the horizon condition (G)/V — 0, where 
the expectation-value is calculated in the ground 
state of Ti. e f / . The horizon condition contradicts 
the usual perturbative expansion, but it is con- 
sistent with an expansion in powers of g$ and the 
Ansatz H 



(M- l ) ab {x,y;A) = g^ 1 S a ' b u(x - y) + 



,(2) 



where M. = —D-V is the 3-dimensional Faddeev- 
Popov operator that appears in TL CO ui ■ To low- 
est order in go , the Schwinger-Dyson equation for 
the Faddeev-Popov propagator becomes an equa- 
tion for the function u which, in the continuum 



limit, reads in momentum space [jjj 
1 n 



u(q) 



(2ir) 



d 3 kD°(k)\ [q 2 -{q- kfk- 



x [u(k) — u(k + q) 



(3) 



where D°(k) is the equal-time gluon propagator, 
to zeroth order in go. 

The kernel D°(k) is an expectation- value cal- 
culated in the ground state ^0 [u] of the effective 
Hamiltonian ~H e f / to zeroth order in g . With ne- 
glect of dynamical fermions, two terms contribute 
to 



a_/0 t_/0 I o_/0 

n eff — n ho "T ri c i 



(4) 



a harmonic-oscillator-like term and a Coulomb- 
interaction one, 



.3 T—l 



d 3 kp a (-k)v(k)p a (k) . (5) 



Here v(k) = k 2 u 2 (k) 7 and p a (k) is the Coulomb- 
gauge color-charge density defined in pl§|. We 
cannot hope to solve this problem exactly be- 
cause it involves the color-Coulomb interaction 
of dynamical gluons. We use instead the ground 
of the harmonic oscillator Hamilto- 



state VP 



(Q)r 




nian only and we obtain 0] 



2w(fc) 2^/k 2 + n*u(k) 



(6) 



"Poster presented by A. Cucchieri. 



An estimate of the leading correction gives a 
renormalization of the length scale, but does not 
otherwise change D°(k) qualitatively. 
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Let w(k) = v,(°\k) be the solution to (^) with 
this kernel. If we define [in the SU (n) case] 



w(q,fj,) 



i 



then (||) becomes 
1 

8^ 
fc 3 (l 



1/2 ; 



q n 



(7) 



<r x (<7) 



rffc 



dz 



with p 2 : 
§ that 

<?(<?) = 

B = 



+ A: 2 



g(fc) _ g(p) 
fc 2 p 2 



(8) 



2kqz. It is easy to check 



£<T 4 / 3 

T(16/3) ^ 2 

r(8/3)r(2/3) 



2/3 



(9) 
(10) 



is a self-consistent solution of the integral equa- 
tion (||) in the IR limit. In jl]|| it has been proven 
that, as q — + oo, the function g(q), defined in (]?]), 
is given by 



l (q) = n(e^)- 1 [log(g 2 mr 2 ) 
+ 3 _1 loglog(g 2 m~ 2 )] 



(11) 



where mi and mi are unknown constants. 

The function g 2 (q) plays the role of a running 
coupling constant 4 7ra(g) and has the UV be- 
havior predicted by the perturbative renormaliza- 
tion group, although the coefficient &o = n/(6ir 2 ) 
does not have the expected value (11 n)/(48 tt 2 ). 
As explained in 0] , one recovers the correct value 
for this coefficient from terms in Tl e ff of higher 
order in go. 

2. A TRIAL SOLUTION FOR THE IN- 
TEGRAL EQUATION 

In order to find an approximate numerical so- 
lution for the integral equation (||) we follow the 
approach in 0] and simplify the problem by using 
a trial solution depending on a set of parameters 
H. We have adopted the form 



?(<?) = 



7T\/6 



{q/qo 



(2+v) 



(lo/q) 



(2+f) 



log(l + g 2 mr 2 ) + S^ 1 log 



+ log (1 + ? 2 m 2 " 2 ) 



-1/2 



(12) 



which has the asymptotic behaviors (||) and (|T] 
This function is positive, monotonically decreas- 
ing and depends on the four parameters mi, 7712, 
go and v. 

We tuned these parameters by evaluating (||) 
numerically using the set / of values of q described 
below. To compare the results obtained for differ- 
ent sets of parameter values we used the quantity 



M{I) = max Ag(q) 
gel 



max 



1 - 



9out{q) 
9in(q) 



(13) 



Here gi n {q) is used on the r.h.s. of (g) and g ut{q) 
is the result obtained on the l.h.s. The set I con- 
sisted of 105 points in the interval [10 _7 ,10 6 ]. 
The accuracy of the numerical integration was 
fixed to five parts in 10 3 . The best values we 
found, q = 1.87, mi = 1.54, m 2 = 2.74 and 
v = -0.06, give M(I) = 0.01826. We also 
evaluated A(I) = (Ag(q)), where the average 
is taken over all the points in /. We obtained 
A(I) = 0.00752 ± 0.00054. 

In order to test our solution we checked [|: (1) 
that the theoretical behavior at small q is sat- 
isfied; (2) that our result is independent of the 
set of points I used for tuning the four parame- 
ters; and (3) that our solution is stable, namely 
if we use different values for the parameters, the 
corresponding output points g ut(q) move in the 
direction of our solution. 

3. THE POTENTIAL AND THE FORCE 

The potential v(k) which appears in (^|) has the 
Fourier transform (apart from an additive con- 
stant) 



V(r) = 



dkg 2 {k) 



sin ( k 1 



n 2 - 1 1 



2n 2 2tt 2 J " a K "' kr 



.(14) 



This gives a force f(r) = — d r V(r) that at 
large separations goes as 



i 2 -l 1 10B 2 r 2 / 3 



2n 2 4tt 3r(ll/3) ' 



(15) 
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and a potential energy that grows as r 5 / 3 . 

From ( [Tl| ) we obtain the limiting behavior of 
the force at small r, 



2 „2\-l 



3w(n- 1) {An z r z ) 



log(A|r 2 ) 
l 



3" 1 log log (A 2 r 2 



, (16) 



where A# = e 7 mi, 7 is the Euler constant, 
and A is a constant. 

4. RESULTS AND CONCLUSIONS 

The potential V(r) which we have calculated 
appears in the quantum field theoretic Hamilto- 
nian Ti e f f . A related quantity is the ground state 
energy E(r) of 7i e // in the presence of a pair 
of external quarks at separation r, H. e ffty = 
E(t)^!q, which is gauge-invariant. They differ by 
the QCD analog of vacuum polarization. How- 
ever because of asymptotic freedom we expect 
that V(r) behaves like E(r) for small r, once a 
physical length scale is adopted. We also expect 
that V{r) and E(r) differ significantly at large 
r. The quantity E(r) is known phenomenolog- 
ically in the range 0.2 fm < r < 0.8 fm. It is 
believed, consistent with lattice gauge theory cal- 
culations [||, that in the absence of dynamical 
gluons E(r) = Kr, holds asymptotically at large 
r, where K is the string tension. The asymptotic 
behavior V(r) ~ r 5 / 3 , is the color-electric field 
energy of two superposed spherically symmetric 
color-electric fields of a quark and anti-quark at 
separation r. Although the power 5/3 may be an 
artifact of the approximation, it is to be expected 
that, for any power that exceeds unity, the ground 
state wave function ^o(A tr ) adjusts itself so that 
at large r the color-electric field is contained in a 
flux tube, which gives a lower energy, that rises 
linearly with r . 

As a test of these ideas and of the approxima- 
tions made, we shall directly compare our results 
for V(r) with phenomenological fits to E(r), to 
see if there is a range of "small" r for which V(r) 
behaves like E(r). To this end, we consider two 
phenomenological models: the Cornell potential 
|?| and the Richardson potential || , which give a 
good fit to the cc and bb spectra. 



From the analytic expression (|l2|), we evaluate 
(numerically) V(r) and fir). To make a connec- 
tion between these dimensionless quantities and 
the real world, we fix the length scale by using 
Sommer's dimensionless phenomenological rela- 
tion H 



rg/(ro) 



-1.65 



(17) 



which holds for the Cornell force at tq a = i?o = 
0.49 fm = 2.48 (GeV)- 1 . 

From the value of a in (GeV) -1 we obtain the 
string tension || 



a = y / min[-/(r)]a 1 



518 MeV 



(18) 



[Because —fir) increases like r 2 / 3 when r goes 
to infinity, we cannot use the standard defini- 
tion a = lmv^oo — fir). Equation ( p"8[ ) defines 
a where / (r) =0, i.e where the potential is ap- 
proximately linear.] It is not easy to estimate an 
uncertainty for the string tension. However, its 
value seems to depend very weakly on the values 
of the parameters of our trial solution &. 

If we identify the parameter A^ in (|16|) with 
the corresponding physical parameter |10|| , then 
from 



A 



MS 



A 



r e 



-7+1-31/66 



(19) 



we obtain A^-§- = 124 ± 12 MeV (see g). 

In Figure Pi we plot our result for f(r), the 
Cornell forces and the Richardson one. Our 
force gets its maximum value for a separation of 
about 2.5(GeV) _1 and it is almost constant up to 
4(GeV) _1 w 0.8 fm, the variation being of order 
12%. 

If, instead of (n7|), we use 



rlfiro) 



= -1.35 



(20) 



which holds for the Richardson force at Rq = 
2.48(GeV)-\ we obtain </W = 468 MeV, A-jm = 
118 ± 12 MeV and the plot shown in Figure [l|b. 
In this case the agreement is even better: in fact, 
our force reaches its maximum value for a sepa- 
ration of about 2.75(GeV)~ 1 , and its variation at 
4(GeV)" 1 is of order 8%. 

The fit is surprisingly good. The force is 
in qualitative agreement with phenomenological 
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models, and the values obtained for the string 
tension and A-p-g are in semi-quantitative agree- 
ment with lattice Monte Carlo and phenomeno- 
logical determinations (see |(|[7j and 111]]). 

It would appear that, with Sommer's normal- 
ization ( |l7| ) or (^0|), the approximate equality 
V(r) « E(r) extends to the range r < 0.8 fm, 
and moreover that the approximations made in 
our calculation of V(r) do not qualitatively de- 
stroy this agreement. Although there is no a pri- 
ori reason to expect that vacuum polarization of 
gluons should not be important in this range, this 
may not be so surprising after all. For although 
vacuum polarization of quarks does "break" the 
string, this is not yet manifest for r < 0.8 fm. 
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Figure 1. Plot of: (i) our force f(R) in the SU{3) 
case (the curve which is decreasing at large R), 
(ii) the forces (for the cc and the bb cases) de- 
rived from the Cornell potential (the two curves 
very close to each other) , and (iii) the Richardson 
force (the highest curve). In the first case we set 
the length scale by using Sommer's proposal (17), 
while in the second case we use the condition ( 2C ) . 



